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We present a scalable method for the tomography of large multiqubit quantum registers. It acquires informa- 
tion about the permutationally invariant part of the density operator, which is a good approximation to the true 
state in many, relevant cases. Our method gives the best measurement strategy to minimize the experimental 
effort as well as the uncertainties of the reconstructed density matrix. We apply our method to the experimental 
tomography of a photonic four-qubit symmetric Dicke state. 
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Because of the the rapid development of quantum experi- 
ments, it is now possible to create highly entangled multiqubit 
states using photons [1-5], trapped ions [6], and cold atoms 
J7|]. So far, the largest implementations that allow for an in- 
dividual readout of the particles involve on the order of 10 
qubits. This number will soon be overcome, for example, by 
using several degrees of freedom within each particle to store 
quantum information [8]. Thus, a new regime will be reached 
in which a complete state tomography is impossible even from 
the point of view of the storage place needed on a classical 
computer. At this point the question arises: Can we still ex- 
tract useful information about the quantum state created? 

In this Letter we propose permutationally invariant (PI) to- 
mography in multiqubit quantum experiments [9]. Concretely, 
instead of the density matrix q, we propose to determine the 
PI part of the density matrix defined as 



Qpi 



1 ^ 



(1) 



where Hk are all the permutations of the qubits. Reconstruct- 
ing gpi has been considered theoretically for spin systems 
(see, e.g., Ref. [10]). Recently it has been pointed out that 
photons in a single mode optical fiber will always be in a PI 
state and that there is only a small set of measurements needed 



for their characterization 111 Ull2ll . 



Here, we develop a provably optimal scheme, which is fea- 
sible for large multiqubit systems: For our method, the mea- 
surement effort increases only quadratically with the size of 
the system. Our approach is further motivated by the fact that 
almost all multipartite experiments are done with PI quantum 
states lll-lj, la] . Thus, the density matrix obtained from PI 
tomography is expected to be close to the one of the experi- 
mentally achieved state. The expectation values of symmetric 
operators, such as some entanglement witnesses, and fidelities 
with respect to symmetric states are the same for both density 
matrices and are thus obtained exactly from PI tomography 
||2|-|4|]. Finally, if qpj is entangled, so is the state g of the sys- 



tem, which makes PI tomography a useful and efficient tool 
for entanglement detection. 

Below, we summarize the four main contributions of this 
Letter. We restrict our attention to the case of N qubits - 
higher-dimensional systems can be treated similarly. 

1 . In most experiments, the qubits can be individually ad- 
dressed whereas nonlocal quantities cannot be measured di- 
rectly. The experimental effort is then characterized by the 
number of local measurement settings needed, where "set- 
ting" refers to the choice of one observable per qubit, and re- 
peated von Neumann measurements in the observables' eigen- 
bases llal . Here, we compute the minimal number of mea- 
surement settings required to recover gp\. 

2. The requirement that the number of settings be minimal 
does not uniquely specify the tomographic protocol. On the 
one hand, there are infinitely many possible choices for the 
local settings that are both minimal and give sufficient infor- 
mation to find qp\. On the other hand, for each given setting, 
there are many ways of estimating the unknown density oper- 
ator from the collected data. We present a systematic method 
to find the optimal scheme through statistical error analysis. 

3. Next, we turn to the important problem of gauging the 
information loss incurred due to restricting attention to the PI 
part of the density matrix. We describe an easy test measure- 
ment that can be used to judge the applicability of PI tomog- 
raphy before it is implemented. 

4. Finally, we demonstrate that these techniques are viable 
in practice by applying them to a photonic experiment observ- 
ing a four-qubit symmetric Dicke state. 

Minimizing the number of settings. We will now present 
our first main result. 

Observation 1. For a system of N qubits, permutationally 
invariant tomography can be performed with 
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local settings. It is not possible to perform such a tomography 



with fewer settings. 

Proof. First, we need to understand the information obtain- 
able from a single measurement setting. We assume that 
for every given setting, the same basis is measured at every 
site 11411 . Measuring a local basis {|0i), |<fe)} is equivalent 
to estimating the expectation value of the traceless operator 
A = 1 0i } (0i | — | ^2} (02 1- Merely by measuring A® N , it is 
possible to obtain all the N expectation values 

<(i4®(*- n > ® l® n )pi), (n = 0,...,N-l) (3) 

and, conversely, that is all the information obtainable about 
qpi from a single setting. 

Next, we will use the fact that any PI density operator can 
be written as a linear combination of the pairwise orthogonal 
operators {X m <g> Y® 1 <g> Z® m <g> l® n ) P i, where X, Y and Z 
are the Pauli matrices. We consider the space spanned by these 
operators for one specific value of n. Simple counting shows 
that its dimension is ZW-n)- The same space is spanned by 
V (N _ n) generic operators of the type (A®^-^ ® l® n ) m . 
We draw two conclusions: First, any setting gives at most one 
expectation value for every such space. Hence the number of 
settings cannot be smaller than the largest dimension, which 
is T>n- Second, a generic choice of T>m settings is sufficient 
to recover the correlations in each of these spaces, and hence 
completely characterizes g-p\. This concludes the proof [15]. 



The proof implies that there are real coefficients 



„(M,' 



such that 



{{x® k <g> y® 1 ® z® m ® t® n ) m ) = 

f^cf' l < m) {{Af N - n) ®t® n ) vl ). (4) 
i=i 

We will refer to the numbers on the left-hand side of Eq. (0|i as 
the elements of the generalized Block vector. The expectation 
values on the right-hand side can be obtained by measuring 
the settings with Aj for j — 1,2,..., T>n- 

Minimizing uncertainties. We now have to determine the 
optimal scheme for PI tomography. To this end, we define our 
measure of statistical uncertainty as the sum of the variances 
of all the Bloch vector elements 



(<£total) — 

53 £ 2 [{X® k ®Y® 1 ® Z® m ® t® n ) m ] 

k+l+m+n=N 

N\ 
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(5) 



where the term with the factorials is the number of different 
permutations of X® k <g> Y® 1 <g> Z® m <g> 1®". Based on Eq. ©, 
the variance of a single Bloch vector element is 



£ 2 [{X m ®Y 
v 



1 ^Z® m ®l® n ) m ) 

£ 2 [(Af N ' n) ® 1®») PI ]. (6) 



= y^ | c ( fe .'. m )i2p2 



Equation (0 can be minimized by changing the Aj matrices 



and the c 



coefficients. We consider the coefficients first. 



Jk,l,: 
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For any Bloch vector element, finding 
mize the variance Eq. © subject to the constraint that equality 
holds in Eq. dU is a least squares problem. It has an analytic 
solution obtained as follows: Write the operator on the left- 
hand side of Eq. (|6]i as a vector v (with respect to some basis). 
Likewise, write the operators on the right-hand side as Vj and 
define a matrix V = [vi,V2, ■ ■ ■ , Vt> n }- Then Eq. (0|l can be 



cast into the form v = Vc, where c is a vector of the c- ' 
values for given (k,l,m). If E is the diagonal matrix with 

2U A®( N ~ 



entries E 2 j 



= £ 2 U] 



1 ")pi], then the optimal so- 



v, where the inverse is 



lution is c = E- 2 V T (VE- 2 V T y l 
taken over the range [16]. 

Equipped with a method for obtaining the optimal c- ' 's 
for every fixed set of observables Aj , it remains to find the best 
settings to measure. Every qubit observable can be defined by 
the measurement directions 3j using Aj — a,j tX X + ctj t yY + 
a,j tZ Z. Thus, the task is to identify T>n measurement direc- 
tions on the Bloch sphere minimizing the variance. In gen- 
eral, finding the globally optimal solution of high-dimensional 
problems is difficult. In our case, however, £totai seems to pe- 
nalize an inhomogeneous distribution of the Sj vectors; thus, 
using evenly distributed vectors as an initial guess, usual min- 
imization procedures can be used to decrease £total an d obtain 
satisfactory results [16]. 



The variance £ 2 [(A 



®(N-n) 



l® n )pi] of the observed 



quantities depends on the physical implementation. In the 
photonic setup below, we assume Poissonian distributed 



counts. It follows that (see also Refs. B171I18IO 
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where (AA) 2 g = (A 2 ) e - (A) 2 g , g is the state of the sys- 
tem and Aj is the parameter of the Poissonian distribution, 
which equals the expected value of the total number of counts 
for the setting j. The variance depends on the unknown state. 
If we have preliminary knowledge of the likely form of go, 
we should use that information in the optimization. Other- 
wise, go can be set to the completely mixed state. For the 
latter, straightforward calculation shows that £ 2 [(A, ® 

t®") PI ] = (^y /(Aj - 1). For another implementation, 
such as trapped ions, our scheme for PI tomography can be 
used after replacing Eq. ^} by a formula giving the variance 
for that implementation. 

Estimating the information loss due to symmetrization. It 
is important to know how close the PI quantum state is to the 
state of the system as PI tomography should serve as an alter- 
native of full state tomography for experiments aiming at the 
preparation of PI states. 

Observation 2. The fidelity between the original state and 
the permutationally invariant state, F(q,Qpi), can be esti- 
mated from below as F(g, gp\) > (P s )g, where P s — 



(a) 



(c) 
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Figure 1: (Color online) (a) Comparison of the 34 symmetrized 
correlations coming from (crosses with error bars) 15 permutation- 
ally invariant measurement settings with optimized Aj matrices for 
N — 4 qubits and (diamonds) from full tomography requiring 81 
local settings. The average uncertainty of all symmetrized correla- 
tions obtained from full tomography is ±0.022, and is not shown in 
the figure. The labels refer to symmetrized correlations of the form 
given in the left-hand side of Eq. l|4}. The results corresponding to the 
15 full four-qubit correlations are left from the vertical dashed line, 
(b) Measurement directions. A point at (a x , a y , a z ) corresponds to 
measuring operator a x X + a y Y + a z Z. (c) Results for randomly 
chosen Aj matrices and (d) corresponding measurement directions. 
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J2 n =o \^n )(Dn \ i s me projector to the TV-qubit symmet- 
ric subspace, and the symmetric Dicke state is defined as 

\ D N l) ) = CP EMM® {N ~ n) ® |l)^ n ), where the 
summation is over all the different permutations of the qubits. 
Observation 2 can be proved based on Ref. [ 19] and elemen- 
tary matrix manipulations. Note that Observation 2 makes it 
possible to estimate F(g, Qpi) based on knowing only gpj. 

Lower bounds on the fidelity to symmetric Dicke states, i.e., 
Tr(\D^'){D^'\g) can efficiently be obtained by measuring 
X, Y and Z on all qubits, i.e., measuring only three local set- 
tings independent of N B20I1 . With the same measurements, 
one can also obtain a lower bound on the overlap between the 
state and the symmetric subspace. For four qubits, this can be 
done based on P s > [(J* + J 4 + J*) - {J* + J 2 y + J z 2 )]/18, 
where J x = (l/2)£ fc X fc , J y = (1/2) J2 k Y k , etc. Opera- 
tors for estimating (P s ) for N — 6, 8 are given in Ref. Ilql . 
This allows one to judge how suitable the quantum state is for 
PI tomography before such a tomography is carried out. 

Experimental results. We demonstrate the method and the 
benefits of our algorithm for PI tomography for a 4-qubit sym- 



metric Dicke state with two excitations \D\ '). First, we opti- 
mize the 3j :'s and the c ■"' 'sforgo = 1/16 and only for the 
uncertainty of full four-qubit correlation terms, which means 
that when computing £totai, we carry out the summation in 
Eq. <(5j only for the terms with n = 0. With simple numerical 
optimization, we were looking for the set of Aj basis matri- 
ces that minimize the uncertainty of the full correlation terms. 
Then, we also looked for the basis matrices that minimize the 
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Figure 2: (Color online) (a) The real and (b) imaginary parts of the 
density matrix coming from full tomography. (c),(d) The same for 
permutationally invariant tomography with optimized and (e),(f) ran- 
dom measurement directions, respectively. 



sum of the squared error of all the Bloch vector elements and 
considered also density matrices different from white noise, 
such as a pure Dicke state mixed with noise. We find that the 
gain in terms of decreasing the uncertainties is negligible in 
our case and that it is sufficient to optimize for go = 1/16 and 
for the full correlation terms. To demonstrate the benefits of 
the optimization of the measurement directions, we also com- 
pare the results with those obtained with randomly distributed 
basis matrices. 

The Dicke state was observed in a photonic system. Es- 
sentially, four photons emitted by the second-order collinear 
type-II spontaneous parametric down-conversion process 
were symmetrically distributed into four spatial modes. Upon 
detection of one photon from each of the outputs, the state 
\D\ ') is observed. Polarization analysis in each mode is 
used to characterize the experimentally observed state. We 
collected data for each setting for 5 minutes, with an average 
count rate of 410 per minute. The experimental setup has been 
described in detail in Refs. Q2L l3|] - 

First, to check the applicability of the PI tomography, we 
apply our tools described above requiring only the measure- 
ment of the three settings, X® 4 , Y m and Z m . We determine 
the expectation value of the projector to the symmetric sub- 
space, yielding (P s ) > 0.905 ± 0.015. Based on Observation 
2, we obtain F(g, g PI ) > 0.819 ± 0.028. These results show 



that the state is close to be PI and has a large overlap with the 
symmetric subspace. Thus, it makes sense to apply PI tomog- 
raphy. 

For PI tomography of a four-qubit system, the measurement 
of 15 settings is needed. We used Eq. © to obtain the Bloch 
vector elements from the experimentally measured quantities. 
This way, we could obtain all the 34 symmetric correlations 
of the form (X® k <g> Y® 1 ® Z® m <g> t® n ) vl . In Fig.Q] we give 
the values of the correlations for optimized and for randomly 
chosen measurement directions, compared to the results ob- 
tained from full tomography, which needed 81 measurement 
settings. As can be seen in Fig. Q] the uncertainty for the op- 
timized settings is considerably smaller than the one for the 
randomly chosen settings. Moreover, the results from the op- 
timized settings fit very well the results of the full tomography. 
In Fig. [2] we compare the density matrices obtained from full 
tomography [Fig. 0a)], from PI tomography for optimized 
[Fig. 0b)] and for random measurement directions [Fig. 0c)]. 
Because of noise, the fidelity of the result of the full tomogra- 
phy with respect to \D\ } is 0.873±0.005, which is similar to 
the fidelity of the results of the PI tomography with optimized 
settings, 0.852 ± 0.009 Q2lb - In contrast, for the method using 
random measurement directions, the fidelity is 0.814 ± 0.059, 
for which the uncertainty is the largest compared to all previ- 
ous fidelity values. Finally, we also computed the fidelity of 
the results with respect to the PI density matrix obtained from 
full tomography [22]. The results of the PI tomography with 
optimized settings shows a good agreement with full tomog- 
raphy, the fidelity is 0.947, which is quite close to the fidelity 
between the results of full tomography and its PI part, 0.964. 
On the other hand, for the PI tomography with random settings 
the corresponding fidelity is much lower, 0.880. Overall, the 
PI tomography shows a good agreement with the full tomog- 
raphy for this particular experiment. However, a reasonable 
choice of measurement directions is needed to obtain uncer- 
tainties in the reconstructed Bloch vector elements similar to 
the ones from full tomography. 

Finally, let us comment on how our method can be extended 
to lager systems. Permutationally invariant operators can be 
represented efficiently on a digital computer in the basis of 
(X® k (Ei Y® 1 <g> Z® m ® l® n ) P i operators. We determined 
the optimal Aj operators for PI tomography for systems with 
N = 6,8, ..., 14 qubits. To have the same maximum uncer- 
tainty of the Bloch vector elements as for the N = 4 case, one 
has to increase the counts per setting by less than 50% 111 611 . 

In summary, we presented a scalable method for permuta- 
tionally invariant tomography, which can be used in place of 
full state tomography in experiments that aim at preparing per- 
mutationally invariant many-qubit states. For our approach, 
the same operator has to be measured on all qubits, which is 
a clear advantage in some experiments. We showed how to 
choose the measurements such that the uncertainty in the re- 
constructed density matrix is the smallest possible. This paves 
the way of characterizing permutationally invariant states of 
many qubits in various physical systems. Moreover, this work 
also shows that, given some knowledge or justifiable assump- 



tions, there is a way to obtain scalable state tomography for 
multiqubit entangled states. 
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Supplementary Material 

The supplement contains some derivations to help to under- 
stand the details of the proofs of the main text. It also contains 
some additional experimental results. 

Proof of that we have to measure the same operator on all 
qubits. From the proof of Observation 1, we know that at least 
T>m measurements are needed to get the expectation values 
of all the £>jv independent symmetric full TV-particle correla- 
tions. What if we measure T>m settings, but several of them 



(1) 4 (2) 
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are not {A,-, A,-, ..., A,-}-type, but {A) 1 ' ,A) Z > , ...,A)" >}- 
type, i.e., we do not measure the same operator on all qubits? 
Each setting makes it possible to get a single operator con- 
taining full A-qubit correlations. Let us denote this operator 
by Mk for k = 1,2,..., T>^. Then, we know the expectation 
value of any operator of the space defined by the Mk oper- 
ators. However, not all M^'s are permutationally invariant. 
Thus, the size of the PI subspace of the space of the Mk opera- 
tors is less than T>n . We do not have T>n linearly independent 
symmetric operators in this space. Thus, T>n measurement 
settings are sufficient to measure g-p\ only if we have settings 
of the type { A,- , A,- , . . . , A,- } . 

Derivation ofEq. (7). The eigen-decomposition of the cor- 
relation term is 



(-4 



®(JV-n) 



1®") PI = J2 Aj,n,fc|*j,fc)(*j,fc|. (SI) 



The individual counts Nc(Aj)k follow a Poissonian distribu- 
tion f(n c , Xj.k), where Aj.fc are the parameters of the Poisso- 
nian distributions and J2k^j-k = ^j- The conditional vari- 
ance, knowing that the total count is Nc{Aj), is 



£: 2 [(Af Ar -" ) ®l^)p I |7V c (A J )] 
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After straightforward algebra, the variance is obtained as 
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(S3) 



Similar results can be obtained through assuming Poisso- 
nian measurement statistics and Gaussian error propagation 
[S1.S2]. If go = l/2 N , then A(Af {N ~ n) ® l® n ) PI is in- 
dependent from the choice of A, . By substituting Aj = Z, 
straightforward calculations gives 



(fc,/,m) 
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(S4) 



Obtaining the formula for c- ' ' for the smallest error. 

We look for c- ' for which the squared uncertainty given in 



Eq. (6) is the smallest. In the following, we use the definition 
given in the main text for c, v, V and E. Thus, V is matrix 
mapping a large space R z to a small space R s . Let E be a 
non-singular diagonal matrix in the small space. We have to 
solve 



\Ec\\ 



s.t. 



Vc 



(S5) 



where | \a\ \ is the Euclidean norm of a. Using Lagrangian mul- 
tipliers, we write down the condition for a minimum fulfilling 
the constraints Vc = v 



?=i 



V s {c*E 2 c + V A, [(Vc), -w l }}= 0. 



(S6) 



Hence, the condition for a local (and, due to convexity, global) 



minimum is 



-E- 2 V T \ 



(S7) 



where A 6 R s is the vector of multipliers. In other words, we 
have a minimum if and only if c G range E~ 2 V T . Because 
the range of V T is an s-dimensional subspace in R , there is 
a unique c in that range such that Vc = v. A solution in a 
closed form can be obtained as 

c = E- 2 V T {VE- 2 V T )- 1 v. (S8) 

Simple calculation shows that the Vc = v condition holds 



Vc = VE-'V 1 [VE-'V 



-2t/T\-1- 



(S9) 



Proof of Observation 2. The eigenstates of J 2 = J 2 + 
J 2 + J 2 are usually labelled by \j, to, a), where J^\j, m, a) — 
j(j + l)\j, to,, a), J z \j, to, a) — m\j, m, a), and a is used to 
label the different eigenstates having the same j and to [S3]. 
Let Pj ia denote the projector to the subspace of a given j and 
a. The number of subspaces is denoted by A^ss, and, for a 
given N, it can be calculated from group theory. Moreover, 
P s = Pjv/2,i- Using this notation, g Pi = Y,j, a P o,<xQ P i,a = 
(P s gP s ) + Y,j<N/2, a ( p j-aQ p j,a)- In the basis of J 2 eigen- 
states, gpi can be written as a block diagonal matrix 



£pi = ® {{ P 0,a) e Qi,a) 



(S10) 



j.« 



where gj tCt are density matrices of size (2j + 1) x (2j + 1). 
In another context, 



gpi - 2 t { p j,ot) e Qj,ot, 



(SH) 



0, a 



where gj >a — Pj ta gPj^ a /Tr(Pj !a gPj^ a ). Based on that, we 
obtain 



F{g,gj, a ) = (Pj, a ) t 



(S12) 



Then, due to the separate concavity of the fidelity, i.e., 

F(g,Pigi + P2£>2) > P\F(g,gi) +P2F(g,g 2 ), we obtain 



Table SI: Fidelities to the 4-qubit Dicke states. 



measurement 


l^ 0) > 


\D?) 


\D?) 


lM 3) > 


l<> 


E 


full tomography 


-0.001 ± 0.002 


0.023 ± 0.004 


0.873 ± 0.005 


0.026 ± 0.004 


0.002 ± 0.002 


0.922 


full tomography (max-like) 


0.001 


0.021 


0.869 


0.023 





0.914 


PI tomography 


-0.001 ± 0.002 


0.040 ± 0.007 


0.852 ±0.009 


0.036 ± 0.007 


-0.002 ± 0.002 


0.925 


PI tomography (max-like) 


0.003 


0.038 


0.850 


0.037 





0.928 


PI tomography (ran) 


0.000 ± 0.002 


0.055 ± 0.027 


0.814 ±0.059 


0.023 ±0.027 


0.001 ± 0.002 


0.893 


PI tomography (ran,max-like) 


0.004 


0.050 


0.816 


0.020 


0.007 


0.897 
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Figure SI: (a) The difference of the real part of the density matri- 
ces from optimized settings and the one of full tomography, (b) The 
difference of the density matrices from random settings and the one 
of full tomography. For the former, no clear structure is observed, 
whereas for the latter the largest difference is observed for the antidi- 
agonal elements. 



F(q,qpi) > 
Substituting 
F(g, gpi) > 

\Ps)g ± 2^j<JV/2 



{P B ) e F(Q,Q B ) + 

Eq. |[ST2t into 



P s/g + l^j<N/2,a\Pj,a/g- 



\Pj,a)t 



}-ij<N/2,a\Pj, a ) e-^ \Q> Qj,oc)- 

this inequality, we obtain 
Using the fact that 
1, we obtain 



F(g,gpi) > (P s )l + ( ' 



<^> e ) 2 



iVss-1 



(SB) 



In many practical situations, the state g is almost symmetric 
and N is large. In such cases the second term in Eq. JS 13b 
is negligible. Thus, a somewhat weaker bound presented in 
Observation 2 can be used. 

Numerical optimization used to minimize £totai- The 
measurement directions minimizing £ total can be obtained as 
follows. Let us represent the measurement directions by three- 
dimensional vectors {<%} J^. The operators can be obtained 
as Aj = aj X X + a,j. y Y ± a,j yZ Z. 

First, we need an initial guess. This can come from a set 
of randomly chosen vectors representing the measurement di- 
rections. One can also use the result of a minimization for 
some measure that characterizes how equally the vectors are 
distributed. Such a measure is defined by 



H{»iY) = X>* • %? 



(S14) 



where {?/. represent the measurement directions and • is the 
scalar product and m is an integer. Such cost functions, called 
frame potentials, appear in the theory of i-designs essentially 
for the same purpose. 

After we obtain the initial guess from such a procedure, 
we start an optimization for decreasing £ total- At each itera- 
tion of the method, we change the measurement directions by 
rotating them with a small random angle around a randomly 
chosen axis. If the change decreases £totai, then we keep the 
new measurement directions, while if it does not then we dis- 
card it. We repeat this procedure until £total does not change 
significantly. 

Three-setting witness for estimating the fidelity The 
three-setting witness for detecting genuine multipartite entan- 
glement in the vicinity of the Dicke state is [S4] 



W 



(P3) _ 
D(4,2) 



= ^ + \( J l+ j2 v- J t-Jt) + TlJl-T-2 J i- ( S15 ) 



For this witness we have [S4] 



W, 



(P3) 



3W 



(P) 



K D(4,2) °'' , 'D(4,2) — "' 

where the projector witness is defined as 



W> 



(P) 

D(4,2) 



i-\d™)(d[ 



(S16) 



(S17) 



,(2) 



Hence, the fidelity with respect to the state \D)£' ) is bounded 
from below as [S4] 



F, 



D(4,2) 



> 



±(W (P3) ) 

3\^ V D(4,2)/ 



(S18) 



Fidelities with respect to the four-qubit Dicke states. In 

Table lSll we summarize the results for full tomography (full) 
and for permutationally invariant tomography (PI) for random 
(ran) and optimized (opt) directions. To obtain a physical 
density matrix with non-negative eigenvalues we perform a 
maximum-likelihood fit (max-like) of the measured data. In 
Fig. [STJ the differences between the density matrix obtained 
from full tomography and the ones obtained from permuta- 
tionally invariant tomography can be seen. 

Efficient representation of permutationally invariant 
operators on a digital computer. Every PI operator O can 
be decomposed as 



k+l+m+n=N 



,y« 



z 



®m 



JPI- 

(S19) 




10 12 14 



Figure S2: The maximum uncertainty of the Bloch vector elements 
defined in Eq. dS2U for the optimal measurement settings as a func- 
tion of the number of qubits, N, for N — 4,6, 8, 10, 12 and 14. 



Such a decomposition for operators of the form (A®( N 
l® n )pi with A = a x X + a y Y + a z Z is given by 



J2 a H a 



k„lm( k + l + m ) 1 



k.l.rn 



k\l\m\ 



(X® k <g> Y® 1 <g> Z® m <g> 1 ")pi, 

(S20) 
where the summation is carried out such that k + l + m + n = 
N. 

Results for larger systems. We determined the optimal Aj 
for PI tomography for N = 4, 6, ..., 14. In Fig. lS2l we plot the 
maximal uncertainty of the Bloch vector elements 

emax = max £[(X® k ® Y® 1 <g> Z® m ® l® n ) P i] (S21) 

for the total count realized in the experiment Xj = X — 2050 
as a function of N, when the state of the system is go = 1/2 . 
It increases slowly with TV. Thus, for large TV the number of 
counts per measurement setting does not have to increase very 
much in order to keep the maximal uncertainty of the Bloch 
vector elements the same as for the N = 4 case. In particular, 
for N = 14, a total count of 2797 per setting yields the same 
maximal uncertainty as we had for the N = 4 case. 

An upper bound on the uncertainty of PI tomography for 
go different from the white noise can be obtained by using 
[A(A® {N ~ n) ®b® n )pi] 2 go = 1 for error calculations. Accord- 
ing to numerics, for optimal Aj for N = 4, 6, ..., 14, e max re- 
mains the same as in the case of white noise, since for the full 
correlation terms with n = the upper bound equals the value 
for white noise, and the full correlations terms contribute to 
the noise of the Bloch vector elements with the largest uncer- 
tainty. Thus, the total count per setting will not increase more 
with the number of qubits even for states different from the 
completely mixed state. 

The operators that give a bound on (P s ) with three settings 
for N = 6 and 8 are the following 

^ (6) > J5(Q 2 + Jf)-^(Q 4 + J5) + 4± 5 (Q 6 + J z 6 ) ; 

P s (8) > -0.001616Q 2 + 0.002200Q 4 -0.0006286Q 6 

+ 0.00004490Q 8 + 0.003265 J\ - 0.004444 J 2 4 

+ 0.001270 J z 6 - 0.00009070 J 8 , (S22) 



where Q n = J" + J". They were determined using semi- 
definite programming, with a method similar to one used for 
obtaining three-setting witnesses in Ref. [S4]. They have an 



(3)\ 



>W\ 



expectation value +1 for the Dicke states [£>g ; ) and \D% }, 
respectively. Moreover, their expectation value give the high- 
est possible lower bound on (P s ) for states of the form 



g noisy (p)=P^ + (l-p)\D^ /2) )(D 



(S23) 



among the operators that are constructed as a linear combi- 
nation of the operators J, n . The validity of the relations in 
Eq. dS22| ) can easily be checked by direct calculation. 

Bounding the differences between elements of g and gpi 
based on the fidelity. For any pure state |\I/), it is possible to 
bound the difference between |(*|g PI |^)| and |(\I/|£|\I>)| as 

|(*|e|*) - (*|gpi|*>| < V~l - F(g, g P1 ). (S24) 

Thus, if the fidelity is close to 1, then (*|g|*) ~ (*|gpi|*), 
even if \^} is non-symmetric. If |\t) is an element of the prod- 
uct basis, e.g., 1 0011 ), then Eq. JS24t is a bound on the dif- 
ference between the corresponding diagonal elements of g and 

QPl- 

Eq. ( IS24I ) can be proved as follows: There is a well-known 
relation between the trace norm and the fidelity [S5] 

^Hc-CPllltr < y/l-F{g,epi). (S25) 

Moreover, for a projector P and density matrices g^ we have 
[S6] 

\Tr(P Ql ) - Tr (Pg 2 )\ <~\\ Bl - g 2 \\tr- (S26) 

Combining Eq. dS25T l and Eq. (|S26] i, leads to Eq. (|S24l . 
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